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Abstract 

Frequency standards based on atomic states, such as Rb or Cs vapors, or single trapped ions, 
are the most precise measures of time. Here we introduce a complementary device based on spins 
in a solid-state system — the nitrogen-vacancy defect in single crystal diamond. We show that 
this system has comparable stability to portable atomic standards and is readily incorporable 
as a chip-scale device. Using a pulsed spin-echo technique, we anticipate an Allan deviation of 
a y = 10 -12 r -1 / 2 with current photoluminescence detection methods and posit exceeding 10 -14 
with improved diamond material processing and nanophotonic engineering. 
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Atomic clocks are the most accurate systems for measuring time and frequency. They 
are used in a broad range of applications, including communication, computing systems 
and navigation, such as the global positioning system (GPS). Modern frequency standards 
derive their stability from the precisely measured internal hyperfine level splittings of atoms 
of Cs, Rb, or H. When an oscillating magnetic field is resonant with the energy difference 
of these internal states a change in population between levels changes the radiofrequency 
or optical absorption. Standard lock-in techniques modulate the driving frequency and 
monitor the absorption as a correction for a tunable active reference oscillator, e.g. a quartz 
crystal, thus stabilizing it to the atomic line [lj. Recent experiments on single trapped 
ions [2j and ensembles of atoms trapped in optical lattices [31 H] have far exceeded the 
international cesium standard, even enabling the observation of general relativity corrections 
within a few meters [5]. Such precision, however, comes at the expense of mobility, as the 
infrastructure for these standards encompass several tens of cubic meters of space. At the 
other extreme, portable standards based on rubidium vapor cells provide excellent stability 
for time scales ranging from 1 s to 10 4 s and find usage in satellites, laboratory equipment, 
and cellular communications. Mobile devices, which typically do not contain their own 
precision standards, can share GPS time signals for maintaining communication standards, 
but when the external lock signal is obstructed a precise local frequency standard with 
minimal drift is necessary to maintain synchronization. 

To the address this need, several groups have miniaturized these atomic standards on-chip 
through the aid of modern microfabrication techniques applied to detectors and lasers [6]-[8] . 
Here, we propose a solid-state alternative based on electron spin states in the negatively 
charged nitrogen vacancy center (NVC) in diamond. The diamond system offers a host of 
advantages. First, single crystal diamond can be grown into a micron-scale, radiation hard 
chip, which makes it portable and well-suited for integration in a semiconductor fabrication 
process [9j. Second, this solid-state system derives its performance as a clock from excep- 
tional spin lifetimes of the NVC [10] and resembles atomic and molecular systems. The 
optical detection of the NVC also increases the signal-to-noise of solid-state standards based 
on inductive detection [llj. Compared to vapor cells, this standard does not suffer from 
doppler or collisional broadening. Although it exhibits an increased homogeneous linewidth 
compared to atomic standards, due to a complex mesoscopic environment [T2J , a higher 
density of defects in solids allows for comparable frequency stability in smaller sensor vol- 
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umes. We estimate a stability, expressed as an Allan deviation, of a y w 2 x 10 -12 at Is of 
averaging for a 0.1 mm 3 diamond sample. Paired with potentially modest operational power 
requirements, this NV clock could result in a new generation of portable solid-state frequency 
standards. A proposed device is shown in Figure [I] and contains: a diamond chip grown with 
a dielectric cavity for lowering the optical pumping power requirements, photoluminescence 
(PL) collection with on-chip Si photodiodes, and planar microwave waveguides for address- 
ing hyperfine transitions, all of which can be integrated into modern device fabrication. 

532 nm excitation 




FIG. 1: Schematic of the diamond frequency standard. A thin (100 /im) diamond chip is surrounded 
by dielectric stacks (Bragg Reflectors) on both sides to create a resonant cavity for 532nm excitation 
in order to reduce the power requirements. On-chip 532nm excitation comes from a doubled 1064nm 
surface emitting laser (not shown). Silicon photodetectors under the diamond collect emission from 
face. Microwaves, which address the NV magnetic sublevels, are applied to the entire sample by a 
planar stripling. 

An atomic clock derives its stability from the large quality factor, Q = v/Av, of the 
probed resonance, with narrow linewidth, Az/, being much smaller than the resonant fre- 
quency v. In the solid-state, we desire: (i) a microwave transition {y ~ GHz) that also 
exhibits a large Q; (ii) a resonance that does not vary due to material or fabrication pro- 
cesses; and (iii) a precise method of measuring population changes. The NVC in diamond 
satisfies these criteria having a ground state spin triplet with long (>lms) spin coherence 
time (thus a small linewidth); a crystal field splitting of the ground state with an intrinsic 
resonance frequency near 2.8GHz that is independent, to lowest order, of applied magnetic 
field; and spin states that are optically polarizable and measurable on single molecule length 
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scales, due to the spin-dependent relaxation of the defect. 

The magneto-optical description of the NVC is well-documented in the literature, thus 
we give only a brief phenomenological summary. Figure 2A shows the relevant spin sublevels 
(0,1,2,3 and S) for the NVC. Optical absorption of green laser light causes for broadband 
PL of the NVC from 637-800 nm. A spin-dependent intersystem crossing between excited 
state triplet (3) to a metastable, dark singlet level (S) changes the integrated PL for the spin 
states |0) and |±1). The deshelving from the singlet occurs primarily to the |0) spin state, 
providing a means to polarize the NVC. Microwave fields resonant between levels & 1, 
perturb the spin populations, and thus the PL response. This response can be measured in 
a steady-state (or continuous wave) response to simultaneous optical and microwave fields 
(e.g. Figure 2B); or, in a pulsed manner, by preparing a state using only microwaves, and 
observing transient PL response. 

In order to quantitatively compare both methods, we construct a simple model for the 
relevant spin dynamics described by Hamiltonians for the lowest and first excited triplet state 
and two metastable singlet states. For the purpose of a frequency standard, we monitor the 
response of only the ground state triplet sublevels to resonant excitation, yielding the ground 
state Hamiltonian[13j: 



Here are the ground state electric dipole moment components along and perpendic- 
ular two the Cs v symmetry axis of the defect. D gs is the ground state crystal field splitting 
(2.870 GHz), /i& is the Bohr magneton, and g is the Lande factor (assumed to be isotropic). 
Sk are spin-1 operators in the k = {x ) y, z} directions. The local electric field vector, induced 
by crystal strain, is a. In the limit of static magnetic and electric fields much smaller than 
D gS) the eigenfunctions are those of the S z operator, as show in Fig. [2j 

A driving field at frequency u) induces electron spin resonance (ESR) transitions between 
|0) and |±1). On resonance (lj « D gs ), the PL decreases and provides a feedback signal as 
to lock uj to D gs . The dynamics are best viewed as a response to a time- varying magnetic 
field Bi = 2b\ cos(27TUjt)x and transforming H gs into an interaction frame defined by the 
operator V = e 27rzujts ^ . Under the rotating wave approximation, we have: 



H^ s = (D g8 + d\\a z )S 2 z + gfi b S • B + 

d±(T X (S X Sy + SyS X ) + dj_(Ty(S* ~ Sj ) 



H' = (D gs + d\\a z - uj)S 2 z + gfi b B z S z + gfibhS, 



X 



(2) 
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FIG. 2: Nitrogen- vacancy center energy levels and resonant response. A. The lowest lying triplet 
( 3 E, 3 A) and singlet ( 1 £ , ) 1 A) orbital states of the NV~ center. /, 7, ft, and A represent the 
absorption, PL, intersystem crossing, and deshelving rates, respectively. Sublevels and 2 are S z 
eigenstates |0), whereas 1 and 3 are |±1), with the degeneracy lifted by small crystal strain or 
applied magnetic field. The model is simplified by approximating both singlet states as a single 
metastable level. B. The steady-state fluorescence emission of the NVC under continuous optical 
and microwave irradiation, detuned from resonance by A. The PL response function, as derived 
from a master equation approach, is approximated by a Lorentzian: F(A) = — ^f^y+dv 2 ) 1 
where C is the modulation depth and Su is the FWHM fH] . 



The relaxation rates of the excited triplet and singlet states, shown in Fig. [2]A, play 
an important role in the optical pumping and spin measurement [15]. We model the total 
magneto-optical response using a master equation approach: 

P = ^IK> P] + E L *P L l - \ L \LkP ~ \pL\L k (3) 

k 

where p is the density operator for the NVC ground, excited triplet, and effective singlet 
states. The jump operators, have magnitudes corresponding to relaxation rates sjr^. 
Conveniently, the solution to the equation yields the total magneto-optical response for both 
continuous and pulsed excitation, allowing for numerically assessment of clock sensitivity 
with changes to electromagnetic fields. 

Under continuous excitation of optical and microwave fields the NVC frequency standard 
closely resembles a two-isotope Rb standard. Spin-dependent PL of the NVC occurs under 
non-resonant absorption of green light of intensity /. Application of a microwave field of 

intensity Q = gpbbi, detuned from resonance by an amount A = D gs — lj, causes a broad, 

2 

phonon-assisted PL: F(I, f2, A) = 7P22 + ~^P33- Here pf 2 and P33 represent the populations 
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of the first excited state spin sublevels in the steady-state, as analytically derived from the 
master equation description of the system [29j . Figure [2]B shows a typical response for F 
for varied detunings, which displays a Lorentzian lineshape. The stability of the clock can 
be derived from the resonance curve by considering the Allan variance: 

where S/N is the signal-to-noise ratio, which depends on the photon shot noise and imperfect 
modulation of the resonance (i.e. C 7^ 1), and r is the averaging time. The intrinsic 
linewidth is limited by the paramagnetic and nuclear spin environments which fluctuate 
during the measurement. The linewidth broadens if the microwave and optical transitions 
are driven near saturation ('power broadening'), but higher pump powers also increase the 
the depth of the dip (C — >> 1). Far below optical saturation, the PL rate is sufficiently small, 
and modulation depth (C) reduced, to cause a decrease in stability per averaging time. 
As the Q is the most important factor for the stability, the optimal condition is attained 
for values which offset the broadening with the reduced signal. Under these conditions 
we estimate a linewidth of 3.6 MHz, an off-resonance fluorescence rate of ^9400 photon/s 
(accounting for finite detector efficiency), and a 17% modulation depth. We thus obtain 
a y (r) = 8.124 x 10" 5 r' 1 ' 2 for a single NVC. 

In short, the NVC lacks sufficient stability when monitoring the PL response continuously. 
The laser excitation must be reduced far below saturation for the optical power broadening 
of the line to reach the homogeneous linewidth. At such low pump powers, the fluorescent 
photon flux is so small that the gains in Q are offset by losses in S/N . As seen in magnetom- 
etry studies of the NVC [E] , a pulsed microwave excitation scheme that monitors transient 
fluorescence behavior can improve the performance over a continuous excitation/detection 
method. 

Alternatively, the NVC crystal field splitting (D gs ) can be monitored in a pulsed fash- 
ion, akin to Ramsey spectroscopy used in atomic clocks. The standard two-pulse Ramsey 
sequence, with time separation T, imprints a phase proportional to the frequency drift, 5u, 
between two of the hyperfine states of the atomic system. The PL response varies sinusoidally 
with 5ujT and can be linearized to provide a passive standard with frequency uncertainty 
set by X -1 . However, if the system does not remain coherent for times X » X c , the 
modulation of the PL response proportional to the phase does not persist. In atomic vapor 
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cells, trapped ions, or atomic fountains, T c > Is due to magnetic shielding and minimizing 
atomic collisions. In the NV system, due to a complex solid-state environment consisting 
of many nuclear spins coupled to a single defect, T c is limited to tens of microseconds in a 
typical sample. Quantum memories and magnetometers of AC fields employ a Hahn echo 
sequence to extend the coherence time to T 2 by removing slowly varying magnetic fields. 
For frequency standards based on m/ = "clock" states, the additional 7r-pulse of the 
echo sequence completely removes the phase accumulation of frequency drift, making it use- 
less for time-keeping. The S=l nature of the NVC, however, contains symmetries within 
the quadrupolar moments that allow for an modified echo sequence to yield a PL signal 
proportional to the drift for T T 2 . 
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FIG. 3: Spin-1 time keeping pulse sequence. The spin-1 ground state |0) state is prepared by 
optical pumping through the singlet state. 45° rotations ensure that the S z term cancels with the 
echo sequence, while the S 2 Z does not, as illustrated in the projected spherical harmonics (inset). 
A transient fluorescence (TF) measurement records the photocurrent for about 300ns timed with 
a pulse of green light. Equipotential surfaces of the density matrix element for the |0) population 
is plotted according to a rigid-body Euler transformation of the spin-1 operators [TT] . Note that 
the S z term is completely refocused by the echo sequence, while the S 2 operator evolves, resulting 
in a different magnitude between the first and sixth states. 

As seen in Eq. ([2]), the term proportional to S 2 Z is sensitive to the frequency drift (5u) : 
whereas the term proportional to S z varies with B. In addition, hyperfine interactions 
with nearby nuclear spins have a secular correction term to the Hamiltonian oc AuS z I z . By 
modifying the flip-angle of the first and last pulses of an echo sequence from tt/2 to 7r/4, as 
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shown in Fig. [3j the phase accumulation due to the S z term cancels while that from the S 2 Z 
terms add over the duration 2T. 

For simplicity, we monitor the eigenkets of S z to describe the clock, assuming perfect 
microwave pulses and no optical illumination. First, optical pumping of the NV prepares 
the initial state |^ ) = \m s = 0). Evolution under our drift echo sequence, U ec h , gives: 

m = U echo \^ ) = ^sin(0)|O) - i C os(0)(|+l) + |-1)) (5) 

where cj) = (D gs - u)T = 8uT. 

Unlike the steady-state case, the NVC is prepared in a state used for time keeping, then 
measured optically, without simultaneous microwave excitation. The transient PL response 
of the NV center can be modeled using projective measurements. The operator M describes 
the spin expectation value for a PL measurement [18] , M = a|0)(0|+fc(|+l)(+l| + | — 1|) , 
where a and b are independent Poisson random variables. The fractional frequency deviation 
varies as the quantum observable, M for our state \i/jf) according to: 

5uj _ 1 (AM) 

\d(M)/duo\ 

where (AM 2 ) = (M 2 ) — (M) 2 is the variance of the operator. For room temperature spin 
readout of the NVC, 2a = 3b. By calculating the moments of M and assuming we accumulate 
the results for a total time r = M'T, we arrive at: 

/ — \ = C 
\ujq/m> D gs VTr 

with t; w 5 due to a combination of imperfect spin readout (i.e. 6^0), inefficient collection 
of photons (i.e a«l), and a small ratio of shelving state to radiative lifetime (A/7). Near 
the onset of information loss T = T 2 ~ 1ms. With D gs = 2870 MHz, this gives a deviation 
of 5u)/ujo = 8.8 x 10~ 9 /a/t for a single NVC. Thus the pulsed NV frequency standard 
dramatically outperforms the CW scheme. A comparison of the two schemes for different 
defect concentrations, along with canonical standards, is shown in Figure [4} 

The stability for a single NV center can be scaled by collecting the fluorescence of an 
ensemble of TV non-interacting NV centers. A central limit approach to S/N will improve 
the stability by a factor of l/y/~N. To make a reasonable estimate of AT, we start with the 
density of pure diamond: 1.74 x 10 23 C atoms/cm 3 . At an NV - defect fraction of 10 -11 (0.01 
part per billion (ppb) ), the density of defects is 1.74 / /1 m 3 , consistent with commercially 
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available samples. Thus, assuming scaling the standard to an ensemble a 1 mm 3 sample 
would yield a fractional frequency deviation of 2 x 10~ 13 /^/r for the pulsed scheme and 
2 x 10~ 9 /y/r for continuous monitoring. A diamond film of 1mm 2 area, 100 um thick, as 
depicted in Fig. [TJ gives (jP ulsed ^ 2 x 10~ 12 /^. Increasing the defect density to 1 ppb 
would decrease all values by a another order or magnitude, but dephasing effects may start 
to arise due to paramagnetic impurities associated with other defect centers [19] . 

Several deleterious factors may factor into the operation of the NVC clock to further limit 
the performance. 

Dephasing. The coherent evolution of the NVC is limited by interactions with the envi- 
ronment, such that the T -1 / 2 behavior cannot be maintained for all times. Unlike an atomic 
vapor, where collisions and stray magnetic fields dominate the process, the NVC exists in a 
solid-state lattice of paramagnetic impurities and spin- 1/2 carbon nuclei. Excellent material 
growth [10] and precision implantation [20] can yield crystals free of paramagnetic impu- 
rities and provide sufficiently long dephasing times (T 2 ). For a single NVC in an ultrapure 
sample, observable spin coherence after an echo sequence of duration 2T is reduced by a 
factor e -( 2T / T 2) n 5 w here exponent n = 3 due to the nuclear spin bath [2T] . In an ensemble 
of NVC, where each center has a different bath, and thus a different T 2 , the net effect is a 
different time constant, T e , with exponent power between n = 2 and n = 1 depending on 
the NVC density [22] . We note that the NVC defect density in the sample is a fraction of 
the total number of substitutional nitrogen impurities given in the sample. Thus, a heavily 
doped sample, while improving the S/N, can reduce the overall coherence time T e . Based 
on the estimates in Ref. [19] . with an NVC density of 10 17 /cm 3 , a marginal modification can 
be expected. Figure [| shows the Allan deviation floor due to such dephasing considerations. 

Scaling imperfections in ensembles. The pulsed NV frequency standard shows a 
formidable Allan deviation, compared to similarly sized atomic standards, by taking ad- 
vantage of a spin echo. Scaling to N' non-interacting defect centers is the straightforward 
way to reduce the deviation by \fW f . but, several deleterious effects may arise. First, the 
local strain fields within the crystal may vary from center to center due to inhomogeneous 
broadening. We can estimate the effects of strain by assuming a large enough number of 
centers such that Ii z fields given by a zero mean Gaussian distribution with variance S 2 . In 
the limit that this inhomogeneous broadening linewidth exceeds the homogeneous broaden- 
ing, which is responsible for the T 2 time of the coherence, we would further limit the sensing 
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time by a new value T 2 *. We estimate a E « 10 6 V/cm, yielding a T 2 * of « 30ns, according 
to the relationship (T 2 *) -2 = d|£ 2 /2. This would reduce the stability by more than two 
orders of magnitude. Modern CVD samples and annealing techniques may be able to limit 
variations of local strain further. In addition, it is known that the resonance frequency 
corresponding to the zero-field splitting of the NV~ center has a temperature dependence 
dA/dT = —74.2(7) kHz/K [23]. This could be due to lattice expansion and an associated 
strain term that shifts the local strain of each NVC, or it could be to a relaxation of the 
NV Cs v symmetry with temperature dependence. Such temperature dependence makes this 
passive standard behave like a crystal oscillator and poses a challenge for maintaining a pre- 
cise operating temperature. It is not clear how other effects common to crystal oscillators, 
e.g. aging, may be affected by strain relief over a period of years. 




Averaging Time (T) 

FIG. 4: Allan Deviation for Atomic and Solid-state standards. Note that the commercial Rb and 
ensemble NV systems follow the same trend. Note that in moving from a CW scheme with a 
single NV center to an ensemble of centers with pulsed excitation and detection, we gain almost 
six decades of improvement. Legend: x: Al ion clock, o: single NVC, +: ensemble NVC (O.Olpbb 
x 1mm 3 ), commercial Rb clock, □: Rb chip-clock, *:surface acoustic wave oscillator, o:ovenized 
quartz reference, single NVC, steady-state, •: ensemble of NVC, steady-state 

Looking forward, the diamond frequency standard promises a fully chip-integratable solid- 
state time-keeping platform whose performance rivals that of modern atomic clocks, both 
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laboratory and chip-scale. Such miniaturized clocks can be easily integrated into scientific 
and consumer electronics and could be used in harsh external environments. Devices bene- 
fiting include those for wireless communication and GPS navigation with improved tolerance 
to jamming [24]. In addition, the clock's center- frequency of 2.87 GHz would enable wide- 
bandwidth date rate in next-generation cellular communications (/ > 40 GHz). At these 
frequencies, phase-sensitive data encoding is limited by the phase noise of the derived fre- 
quency source, which is usually a MHz range mechanical oscillator frequency upconverted 
by a factor N > 1000. The low noise spectrum of the mechanical oscillator, scales with TV 2 , 
thus degrading performance. A low-power, portable, and stable on-board oscillator in the 
GHz range can avoid this quadratic stability loss (Fig. [!]). 

CONCLUSIONS 

In summary, we have proposed a solid-state spin-optical frequency standard based on the 
NV defect center in diamond. Due to the center's relatively long lifetime, high density of 
spins, and optical detection, we estimate a time stability exceeding a y = 2 x 10 -12 t -1 / 2 , 
rivaling that of the newest chip-scale Cs and Rb standards. Many of the technological 
advances, such as the surface emitting lasers and chip-scale detectors, are readily deployable 
to this system, with the added benefit of the atomic system being inside the diamond 
substrate. We anticipate the ability to improve the frequency standard using repetitive 
readout [25j or entanglement assisted [26] schemes, by using IR absorption detection [27] . 
and by deploying nanophotonic [28] advances in diamond to further improve the sensitivity 
of this device. 
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